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ABSTRACT 

Recent models of molecular cloud formation and evolution suggest that such clouds are dynamic and 
generally exhibit gravitational collapse. We present a simple analytic model of global collapse onto a 
filament and compare this with our numerical simulations of the flow-driven formation of an isolated 
molecular cloud to illustrate the supersonic motions and infall ram pressures expected in models of 
gravity-driven cloud evolution. We apply our results to observations of the Pipe Nebula, an especially 
suitable object for our purposes as its low star formation activity implies insignifcant perturbations 
from stellar feedback. We show that our collapsing cloud model can explain the magnitude of the 
velocity dispersions seen in the 13 CO filamentary structure by Onishi et al. and the ram pressures 
required by Lada et al. to confine the lower-mass cores in the Pipe nebula. We further conjecture that 
higher-resolution simulations will show small velocity dispersions in the densest core gas, as observed, 
but which are infall motions and not supporting turbulence. Our results point out the inevitability of 
ram pressures as boundary conditions for molecular cloud filaments, and the possibility that especially 
lower-mass cores still can be accreting mass at significant rates, as suggested by observations. 
Subject headings: turbulence — methods:numerical — ISMxlouds — ISM:kinematics and dynamics — 
stars:formation 



1. INTRODUCTION 

The supersonic, "turbulent" motions observed in 
molecular clouds must play an important role in star 
formation. Early numerical models of molecular clouds 
often imposed s upersonic velocities as either a continu- 
ous forcing term dMac Low et al1ll998tlStone et al.lll99"l 
iPadoan fc Nord lund 19991) and/or as initial conditions 
(jBate et al.ll2002l 120031 ). Yet without an understanding 
of how supersonic turbulence originates, it is difficult to 
develop a predictive theory of the processes leading to 
the formation of stars. 

The recognition that molecular clouds might of- 
ten if not generally result from accumulation of g as 
by large scale flows ([Ballesteros-Paredes et al.l fl999al fbl: 
lHartmann et al.l I2001D and in p ar ticular from atomic 
flows (IKovama fe Inutsukal 120021 : lAudit fe Hennebelld 
120051 : iHeitsch et al.l l2006bh has made it plausible that 
the turbule nce is arising as a consequence of the cloud's 
formation. IVazquez-Semadeni et al.1 (|2007t ) showed that 
turbulence could develop in clouds formed by variable- 
velocity flows. By imposing a fixed spatial variation in 
the supersonic inflow velocities, they in effect identified 
the d riving mechanism as va riations in the inflow speeds 
(also [H cnnebe lle"et al.| [2008) . In contrast. IHeitsch et all 
( 2006b, 2008b) showed that even uniform inflows can pro- 
duce turbulent substructure if the shock interface is not 
planar and/or precisely perpendicular to the flows. 

It is clear from the studies described in the previous 
paragraph that molecular clouds swept up by supersonic 
flows in the interstellar medium will begin their existence 
as both structured and turbulent. However, this sweep- 
up by itself only results in cold clouds that are at most 
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mildl y supersonic if not subsoni c jKoyama & In utsukal 
| 2002t lAudit fe Hennebeliel 120051 : IHeitsch et al.l l2006bl . 
2008b). The dominant mechanism for producing super- 
sonic "turbulent" motions in molecular clouds - espe- 
cially at colu mn densities typical of star- forming clouds 
- is gravity (lyazquez-Semadeni et al.ll2007l 120081 120091 : 
IHeitsch et a l I l2008bb IHeitsch fe Hartmann! l2008h. This 
was already seen by iBurkert fe Hartmann! ([2004) , who 
pointed out that clouds with many Jeans masses and 
non-spherical geometry are generically susceptible to 
generating large, spatially- variable gravitationally-driven 
flows, as commonly seen in simul ations with non-p e riodic 
gravity allowi n g globa l collapse dBate et al.ll2002l 



gravity allowi n g globa l collapse 1 U ate et al.l zUU/ 
Jta. | Heitsch et all l2008ri IV azquez-Semadeni et al.l 
~7i 120091 ). lHartmann fc Burkertj (|2007l ) went furthei 
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2009D . lHartmann fc Burkertj (|2007l) went further to ar- 
gue that large-scale gravitational collapse is a feature of 
at lea st the Orion A molecular cloud (see also lTobin et al.l 
l2009lL 

Relatively quiescent regions, undisturbed by energy 
input from young stars, can provide good tests of the 
gravity-driven picture of cloud evolution. The Pipe Neb- 
ula is a prominent and well-studied example of a cloud 
of significant mole cular gas mass with little perturb- 
ing st ar formation (lLombardi et al.ll2006l : iForbrich et al.l 
2009). The densest regions of the Pipe generally lie along 
a well-defined filament ary structure (lAlves "etaTJ [20071 
I Rathborne et all [20071: iMuench et al.l 120071 : lLada et all 
120081: IRathborne et al I I2009D - a common feature of 
star-forming clouds (e.g.. ISchneider fc Elmegreenlll979( ). 
While much attention has been paid by the above authors 
to the mass function of the dense cloud cores found along 
the filaments, our interest is in the dynamical environ- 
ment and evo l ution of these cores. 

lLada et "all (|2008l ) concluded that significant external 
pressures were needed to confine many of Pipe cores, 
and attributed these pressures to "the weight of the sur- 
rounding molecular cloud" . As the surrounding cloud 
is highly unlikely to be in hydrostatic equilibrium, this 
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"weight" is most likely to be a dynamic pressure. In par- 
ticular, as the filament and cores represent a significant 
mass concentration, gravitational acceleration of exter- 
nal material should provide a confining pressure. A sig- 
nature of this confining infall should be supersonic line 
broadening. While the internal non-thermal velocity dis- 
persions of the Pipe cor es are subsonic, as detected in - 
dense gas tracers (see also lMversllI983l:lAndre et aL 120071 ; 
iKirk et all [20071 for other regions'). I Onishi et all (jl999l ) 
found supersonic line widths in 13 CO, qualitatively what 
would be expected from gravitationally-accelerated, in- 
falling material which enters the cores and filaments 
through shocks at the boundaries (jGomez et al.1 120071 : 
iGong fe Ostriker| [2009). Our goal here is to make this 
qualitative picture more quantitative. 

In this paper we attempt to develop an understanding 
of the gas flows driven by global gravitational collapse 
onto a dense filament. We start out with an analyti- 
cal model to estimate the expected velocities in the col- 
lapsing gas ($2j). This model is in turn motivated by 
numerical simulations (jj3J| of the flow-driven formation 
of an isolated low-mass cloud. These simulations show 
that filament formation and global collapse are a natural 
consequence of the cloud formation process. We com- 
pare the analytical and numerical results, and then use 
them as a guide to interpret recent observations of dense 
gas filaments in the Pipe Nebula (Sj4]). We find that the 
typical velocity dispersions resulting from collapse of gas 
onto the dense filament agree reasonably well with ob- 
servations of gas at the densities which we probe in the 
simulation. We also find that the ram pressures in the 
numerical simulation agree r easonably w ell with the con- 
fining pressures inferred by lLada et all (2008) for their 
lower-density cloud cores. Our models imply that the 
observed non-thermal velocity dispersions are likely due 
to infall rather than supporting turbulence, and that the 
lower-mass cores can be gaining mass at significant rates. 
We conclude that the gravity-driven picture of cloud evo- 
lution is in reasonable agreement with the observations 
of this quiescent molecular cloud region. 

2. THE INFINITE CYLINDER AND INFALL 

Before examining the results from the numerical sim- 
ulations, it is useful to develop estimates from a simple 
model to serve as a benchmark. Consider a uniform, infi- 
nite, self-gravitating filament extended in the z direction, 
with R the radial distance in cylindrical coordinates. The 
critical line density m c (mass per unit length in the z 
direction), obtained by integrating momentum balance 
from R = to R — oo is a function only of temperature 
()Ostrikerlll96l : 

m c = 2c 2 JG, (1) 

where c s is the isothermal sound speed. Assuming a 
mean molecular weight of 2.36m#, 

m c = 16.3 Tio Mq pc _1 , (2) 

where Tio is the gas temperature in (typical molecular 
cloud) units of 10 K. The filament has a density structure 
as a function of cylindrical radius of 

p = po (1 + i?7(4# 2 ))- 2 , (3) 
where the scale height H is given by 

H = C 2/(2GE ) = 0.19T 1o Vp c - ( 4 ) 



Here we have assumed that the relation between extinc- 
tion and molecular hydrogen column density is Ay — 
1.1 x 10 21 /iV(iJ2). 

We may also relate the scale height to the central den- 
sity; 



H = 



2nGp 



(5) 



For reference, the half-mass radius is at R — 2H . 

Now consider a parcel of gas in free-fall toward this 
cylinder, starting from rest at a cylindrical radial dis- 
tance R - The gravitational acceleration is 

a G = -2Gm/R, (6) 

and thus the infall velocity v at R is 



v = 2(Gmln(R /R)) 



1/2 



(7) 



If we put this in terms of a static filament, 

v^2 3 / 2 c s (m/m c )^ 2 )(\n(R /R))^ 2 (8) 

= 0.53 ln(R /R)) 1/2 {m/m c ) 1/2 T^ 2 kms" 1 . 

The infall velocity is therefore relatively insensitive to 
the initial position. 

In a realistic situation, infalling material will shock at 
the filament boundary as it adds mass to the filament. 
The subsequent velocity dispersion in the p ost-shock, fil- 
amentary gas will be subsonic but non-zero (jGomez et al.l 
[20071: IGong fc Ostriken l2009h . with a precise value de- 
pending upon the rate of radiative cooling. 

This model ignores any gl obal motions which are gener- 
ally present in finite clouds. iBurkert fc Hartmannl (|2004f ) 
showed that the collapse of a circular, uniform sheet re- 
sults in a strong pile-up of material at the infalling edge 
("gravitational edge focusing") 4 . Interior to this edge, 
the collapse timescale t c of a region of extent 5r is rel- 
atively independent of radial position r, and is approxi- 
mately 



tc 



R 



1/2 



(9) 



where R is the inital radius of the sheet and £ is the sur- 
face density. This is also approximately the time taken 
for the edge of the sheet to reach the center. Equation 
© implies that the velocity difference across a region Sr 
is roughly 

*«~(-R-J ( 10 ) 

The relative importance of the global collapse to the 
filament-induced velocities will then depend upon rela- 
tive magnitudes of the surface density of the external 
region and the mass line density of the filament. 

3. NUMERICAL RESULTS 

With the analytical estimates of the flow dynamics 
around filaments in hand, we apply them now to a less 
idealized geometry, using a simulation of cloud growth 
driven by large-scale gas flows. The goal is to test 

4 We s hould point out that this effect has also been noted by 
[Lil J200U1 . who used it to explain clustered star formation, albeit 
in the context of magnetically dominated finite sheets subject to 
ambipolar drift. 
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whether the analytical estimate can reproduce the qual- 
itative behavior of gas collapsing onto a dense filament 
in a more complex environment. We briefly summarize 
the main properties of the si mulation (for further de- 
tails, see iHeitsch et al.l l2008bl ). The simulation models 
the flow-driven formation of an isolated molecular cloud, 
including the appropriate heating and cooling processes 
as well as self-gravity. Two gas flows colliding head-on at 
a shocked interface lead to compression, strong cooling 
and rapid fragmentation due to a combination of dynam- 
ical and thermal effects. The isolated cloud forming out 
of this collision is mildly turbulent, and with increas- 
ing mass it submits to global gravitational collapse per- 
pendicularly to the inflows. The finite cloud geometry 
leads to a sweep-up of material due to global gra vita- 
tional edge focusing (see Burkcrt fe Hartm ann 2004|), re- 
sulting in the forma tion of a dense filament at t he cloud 
edge (see Fig. 1 of IHeitsch fc Hartmann! [200^ ). As we 
will show below, the gas infall onto the filament itself 
is roughly cylindrical. It is the infall of gas onto this 
filament which we are interested in. 

This formation mechanism is able to reproduce some 
of the salient properties of molecular clouds, namely 
their internal turbulence, the predominantly filamentary 
structure of their dense gas, and the observed rapid on- 
set of "star" formation in the clouds. The setup is rather 
generic and is physically equivalent to e.g. the collision 
of two supernova s hells, the sweep- up of gas by an ex- 
panding shell (e.g. JPatel et al. 1998), o r gas swept up in 
spira l arms of galaxies lElmegreen T9791 [20071 : iKim et all 
2003; iDobbs fc Bonnellll2007l) 

With a box size of 22 x 44 2 pc and a resolution of 
256 x 512 2 cells, our simulation does not have the spa- 
tial resolution to directly model the small cores seen in 
the Pipe, since it does not reach the required densities 
and temperatures. Yet, we can use it to check the ap- 
proximate validity of the simplistic model of infall onto 
a filament described in the previous section under more 
realistic conditions, including the cloud's evolution and 
global gravitational collapse. 

To compare the gas infall in our simulation to the 
physically appropriate tracer, we estimate the 13 CO line 
emission from the model cloud. For this estimate, we 
make a simple approximation, motivated by the notion 
that CO formation requires shielding by dust grains (see 
IHeitsch fc Hartmann! [2008 for more details). We decide 
whether CO is "present" in a particular grid cell by de- 
termining the attenuation of the ambient radiation field 
integrated over solid angles. If the effective extinction is 
equivalent to that of an angle-averaged Ay — 1 and the 
local temperature is T < 50 K, we assume CO is present 
in high abundance. Because CO is rapidly diss ociated 
at lower extinctions (Ivan Dishoeck fc Blacki ri988) , we do 
not advect CO for simplicity. The radiation field at each 
grid point is calculated by measuring the incident radi- 
ation for a given number of rays and averaging over the 
resulting sky. The ray number is determined such that 
at a radius corresponding to n c — 256 cells (i.e. half 
the size of the larger box dimensions), each resolution 
element of the cartesian model grid is hit by one ray, 
i.e. n ray ~ 47m 2 .. Thus, fine struc tures and strong den - 
sity variations are resolved (see also lHeitsch et al1 l2006a). 
With the distribution of CO thus determined, we ap- 
ply a three-dimensional version of the radiative transfer 



Montecarlo code by iBernei (|1979f) . as implemented by 
Mardones et al. (in preparation). The code determines 
the level populations by assuming statistical equilibrium 
between collisions and radiation. The collision rates are 
taken from the Leiden Atomic and Molecular Database 5 . 
To lower the computational needs, we reduce the resolu- 
tion of the original data cube by a factor of 2. After few 
tens of iterations, the level populations converge, and the 
line profiles in an arbitrary direction can be obtained. 

Figure [T] shows the velocity-integrated emission 
(grayscale) of the 13 CO-map for a model cl oud assembled 
by larg e-scale colliding flows (model Gf2 of lHeitsch et alJ 
l2008bi at t = 14.5 Myr). The view is along the inflow (x) 
direction, and the frames measure 12 pc across. Over- 
plotted are the spectra for regions with Tf, > 2.5 K. Each 
spectrum covers a velocity range from —2 to +2kms _1 . 
Two features in the line profiles are noteworthy: they are 
frequently asymetric, and some cases they exhibit more 
than a single peak. 

Table[T]lists the first moments of the line profiles. Col- 
umn 1 shows the number of the region, as shown in Fig- 
ure [TJ and column 2 shows the intensity in K kms -1 . 
Columns 3 & 4 list the centroid velocity and the veloc- 
ity dispersion, respectively, and finally, in column 5 we 
report whether the line profile has multiple peaks or not. 

Figure [2] (left) shows the centroid velocity and the ve- 
locity dispersion (right), again in the y — z projection. 
Despite the fact that we view the cloud along the inflow 
direction, the dense filament is coherent in velocity space 
(see also centroid velocities in center column of Tabled]), 
except for the largish core on its right end, which is ap- 
proaching the observer at « —0.5 km s^ 1 . Thus, the ve- 
locity dispersion in the dense gas - the core-to-core veloc- 
ity dispersion - does not contain much information about 
the velocities of the assembling flows (~ 7.9 kms -1 ), in 
fact, the motions in the (coherent) dense post-shock gas 
are subsonic. 

This also can be seen when considering the gas mo- 
tions perpendicular to the filament and in the plane of 
the sky, i.e. the infall of more diffuse gas onto the fila- 
ment (Fig. [3]). The velocity profiles were calculated by 
identifying the filament in the two-dimensional projec- 
tion (see Fig. Q]) . Comparing the velocity and the den- 
sity profiles (top and center panel), it is clear that gas at 
higher densities should show smaller velocity dispersions: 
at a distance of 2 pc from the axis of the filament, the in- 
fall velocity ranges between 0.2 and 0.6 km s" 1 , i.e from 
subsonic to slightly supersonic. In other words, the dense 
filament is essentially a post-shock region (with subsonic 
internal motions), confined by a dynamical pressure gen- 
erated by the (supersonic) infall of gas onto the filament. 
We note that the dense filament is not very well resolved 
and that thus we decided not to model higher-density 
tracers showing the subsonic motions more clearly - 
the line-profiles would not be reliable. Our choice of 
13 CO probes the immediate vicinity of dense cores, which 
makes it an appropriate tracer for our purposes. 

The formation of subsonic cores in post-shock gas has 
also been discussed in the context of driven, supersonic 
turbulence assembling dense co res (jKlessen et alJ l2005t 
iBallesteros-Paredes et al.l l2007f ). In that case, density 
peaks form at locations of maximum compression and of 

5 http: //www. strw. leidenuniv .nl/~moldata/ 
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Fig. 1. — Total intensity map (in K km s 1 ) of 13 CO emission from the model cloud, seen along the inflows. Line-of-sight velocity 
spectra are overplotted for regions with I > 2.5 K km s . The numbers refer to Table[T] 



minimum relative velocity difference. In our case (Fig.[3|), 
the inflows are driven by the deepening gravitational po- 
tential of the assembling cloud, not relying on externally 
forced turbulence. 

The bottom panel of Figure [3] compares the median 
value of the pressures as an estimate of what processes 
dominate the gas dynamics in the environment of the 
filament. The thick solid line corresponds to the ram 
pressure obtained with the infall velocity along the z- 
direction (i.e. perpendicular to the filament in Figure[T]), 
while the thin solid line shows the total ram pressure 
profile. The gravitational pressure p$ is shown by the 
dot-dashed line, while the thermal pressure is given by 
the dashed line. Clearly, the region is more dynamical 
than just indicated by the mean infall velocities: as sug- 
gested by this figure, the infall (upper panel) is driven by 
the gravitational potential (lower panel) , while the ther- 
mal pressure responds only weakly to the gravitational 
pressure. 

To get a clearer view of the dynamics, we show in Fig. 3] 
the three projections of the datacube, with the line pro- 
files overplotted, as in Fig. Q] The x — y panel is on the 



top of the cube, the x — z on the left, and the y — z on 
the right. These projections reveal that the filament is 
actually tilted by approximately 45 degrees in the x — y 
plane, with the left side closer to the observer located 
at x — > oo. From the three-dimensonal datacube, we 
know that the infall velocities in the three directions are 
of similar magnitude. However, this does not necessar- 
ily translates to an infall, double-peaked line profile. Al- 
though the three projections exhibit some double-peaked 
profiles, most of our 13 CO lines exhibit only asymmetric, 
supersonic (a v ~ 0.6 km s -1 ) profiles. Thus, the super- 
sonic line widths observed in molecular clouds that are 
usually attributed to random motions (or " turbulence" ) , 
very likely contain a substantial component of ordered 
motions due to coordinated global collapse. A filament 
tilted with respect to the line-of-sight - a very likely sce- 
nario -, only helps to suppress clear signatures of coor- 
dinated infall onto the filament, giving the impression of 
"turbulent" motions. 

Table [5] summarizes the measured properties of the 
model filament, as identified in the two-dimensional pro- 
jection. Since the massive cores potentially could affect 
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Fig. 2. — Left: Centroid velocity (in km s 1 ) of the model cloud, seen along the inflows. The filament is coherent in velocity space. Right: 
Velocity dispersion (in kms -1 ) of the model cloud, as seen along the inflows. The densest regions have the smallest velocity dispersion. 

the geometry, we show numbers for the full filament (in- 
cluding the cores), and the left half of the filament, which 
does not contain massive cores. 

We note that the filament line density is approximately 
twice the expected limit for the isothermal case (for the 
higher temperature of 28 K in the simulation) . As men- 
tioned above, the filament is tilted, resulting in a shorter 
projected length, and thus in a higher line density. Also, 
other effects such as limited resolution and the possibil- 
ity that the filament is not static but is still accreting 
mass on timescales of order 1 Myr, could play a role. 

Given that the model filament does not agree with the 
static approximation, we use equation Q for the full 




\ What 
2 pc 



Fig. 3. — Velocity, density and pressure profiles of gas falling 
onto the filament, in the plane of sky, against distance to the fila- 
ment. Negative distances are located below the filament. For the 
velocities, the solid line denotes the centroid (mean) velocity, and 
the dashed lines indicate the la scatter. For the density, solid lines 
stand for median values, while dashed lines indicate the lower and 
upper quartiles. The pressure plot shows the ram, internal and 
gravitational pressure as indicated. The thin solid line stands for 
the total ram pressure, i.e. all three components of the velocity 
contribute. 



1 /2 

filament to estimate v ~ 0Aln(R o /R) kms 
value of R to use is not clear, but if we take R a 
and R = 0.2 pc for typical values, v ~ 0.6 kms , this 
is about twice the average infall velocity at the filament 
(Figure EJ. 

To summarize, Figure [3] demonstrates the overall col- 
lapse of the cloud in addition to acceleration onto the 
filament. Obviously, equation (fl~0|) can describe the gas 
dynamics of the simulated cloud only in a very general 
way, as (a) matter is being added to the cloud during 
its collapse from the inflows, (b) the cloud (or filament) 
cross-section is not exactly circular, and as (c) there are 
shear components in the flows. Nevertheless, the overall 
collapse motions are approximately cylindrical (sec total 
ram pressure profile in Fig. [3|), indicating that at this 
stage of the simulation the massive filament dominates 
the gravitational potential and thus the dynamics of the 
cloud, and justifying equation (fT0|) as a rough description 
of the actual dynamics. The total mass of the cloud, in- 
cluding the filament itself, at this epoch is ~ 2000M Q 
distributed over a region approximately 10 x 6 pc, which 
implies an average surface density of £ ~ 7x 10 -3 g cm -2 
or 1.8 x 10 21 cm" 2 . Assuming the initial radius in the 
y-direction was ~ 10 pc, equation (|10p then implies a ve- 
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Fig. 4. — The three intensity projections of the datacube. The flows are along the x-axis, and Figure [T] shows the (y, z)-plane. 



locity gradient across a region of 2 pc of 5v ~ 0.4 km s -1 , 
which again is in reasonable agreement with the results 
in Figure 2. Considering the very crude nature of these 
estimates, agreement with the numerical simulation at 
the factor of two level is adequate. 

4. COMPARISON WITH THE PIPE NEBULA 

Motivated by the global, approximately cylindrical in- 
fall observed in the simulation (@ under rather general 
conditions, we first discuss in this section our analytical 
estimates (|J5J| and then the simulation results in the con- 
text of the Pipe nebula. We summarize the relevant ob- 
servational findings first, and then use the approximate 



agreement between observations and models to interpret 
the physical conditions of the Pipe nebula in the context 
of global infall (fgj). 

Th e most recent es t imate of core properties comes 
from iRathborne et all (|2009h . The core masses range 
from ~ O.3M to about 25M Q , with radii ranging from 
~ 0.04 to ~ 0.2 pc and a typical internal density (H2) of 
- 10 4 cm" 3 . The FWHM of the C 18 lines ranges from 
about 0.4 to about lkrns" 1 . The estimated visual ex- 
tinction rises slowly with mass until about M ~ 2M Q , at 
which point Ay rises rapidly from about ~ 4 to 10 — 15; 
this is roughly the mass range which becomes gravita- 
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number / [K km s 1 ] (v) [km s 1 ] crv [km s ] mult 



TABLE 1 

Parameters of subfields in Figure [T] The last column 

INDICATES THE EXISTENCE OF AN OBVIOUS SECOND COMPONENT. 



M T L m c m c 

[Mq] [K] [pc] measured expected 

full filament 1.1 X 10 3 28 12 9~7) 4~5 

without cores 3.5 X 10 2 27 5.7 6.2 4.4 



TABLE 2 

Filament parameters for models, to be compared to the 

toy model of 53 The mass per length m c is given in 
Mq PC -1 . The first line gives the numbers for the full 
filament (as used for Figure O, while the second line 
refers to the left half of the filament which does not 
contain any cores. 

tionally bound, according to lLada et alJ (|2008h . Median 
properties of the cores are M ~ 0.8M Q , R ~ 0.07 pc, 
and FWHM ~ 0.36 km s -1 . Such a " median" core is typ - 
ically pressure-confined according to lLada et "all {2008) . 
To simplify the comparison with observations, we assume 
a typical temperature of 10 K, though there is some ev- 
idence that some of the Pi pe cores have slightly h igher 
temperatures, - 12 - 15 K (jRathborne et alJl2007t l. 

To proceed further we need an estimate of the average 
line density of the filament. The Pipe, though highly 
elongated, is not a perfectly straight filament nor is it 
unifor m. The sum of the core masses in lRathborne et aTl 
(2009) is 228M Q ; distributing this mass over a length 
~ 9° ~ 20 pc results in an average line density ~ 
llM pc _1 . This estimate neglects the inter-core mass 
but includes the "bowl" region which is much more com- 
plex in structure. If we restrict attention to the por- 
tion of the pipe at negative galactic longitudes, the total 
core mass is 105M Q over ~ 15 pc or m ~ 7M pc~ 1 . 
As the Pipe does not exhi bit considerable star f orma- 
tion at the present epoch (|Forbrich et al.l (|2009|) esti- 
mate a star formation efficiency of 0.06% when compar - 
ing to the total cloud mass given bv lOnishi eta l. 1999), 
it seems reasonable to assume an average line density 
of m ~ 15M pc _1 , to include mass outside the core 
but not so much that the filament would be radially 
gravitationally- unstable. 

To further fix ideas we take an average filament optical 
depth to be Ay ~ 2. This is near the low end of th e 
core extinction values found in iRathborne et al.l {2009) . 
Then from equation (|4]), the scale height is H ~ 0.095 pc; 
this seem s reasonable, seeing that the median core radius 
given by IRathborne et al.l (|2009l ) is about 0.07 pc. We 
may also then derive a central density from equation ([5]) 
of po ~ 9.5 x 10~ 21 gcm~ 3 or n(H 2 )o ~ 2.4 x 10 3 cm~ 3 . 
This density is a few times smaller than the average core 
density n{H2) ~ 10 4 cm~ 3 , but we have not taken any 
confining pressure into account. 

We are now in a position to calculate an estimate of 
the infall velocity using equation (|9|). Taking R /R = 10, 
we find v ~ 0.77 km s . The timescale for this motion is 
sensitive to the choice of initial condition, i.e. the initial 
velocity; for the above solution, the time taken to fall in 
from 0.9 to 0.1 pc is about 2 Myr. 

The translation of this infall velocity into a typical 
velocity dispersion is sensitive to geometry at the fac- 
tor of two level. Assuming that the filament is oriented 
perpendicularly to our line of sight, and the emission 
is simply mass weighted, the average velocity width is 
~ (4/7r) x 0.77kms _1 . If we identify this with the 
FWHM or line width, our results are consistent with 
the median 13 CO line width of ~ 1.2kms _1 found by 
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lOnishi et al.l (|1999l see their Fig. 6b). It is larger than 
that observed in the C 18 cores by about a fa ctor of 2.5; 
however, as pointed out bv lLada et alj (|2008l ). the C 18 
observations sample material both outside and inside the 
cores. As mentioned in §§2 and 3, cores would be pre- 
dicted to be the sites of post-shock gas and thus with 
much smaller internal velocity dispersions. Lada et al 
state that, in the lower-mass cores, half or more of the to- 
tal line-of-sight col umn density arises from material out- 
side the core (also iGomez et al.1 l2007t iGong fe Ostrikerl 
2009). Conversely, this means that for many cores, half 
or less of the emitting material is outside the core. Thus 
we would predict a core velocity width perhaps half that 
above, or aro und ~ 0.4 km s" 1 , in rough agreement with 
observations (|Muench et al.ll2~007f) . In addition, there is 
no specific reason to assume that the filament is oriented 
perpendicularly to the line of sight. Any (very likely) tilt 
of the filament with respect to the line of sight will result 
in an even smaller projected velocity dispersion. 

Going back to the numerical simulations, we note that 
the velocity dispersion in the flow direction, calculated 
as the second moment of the line profile (see Fig. Q] and 
Table [TJ, are or the order of ~ 0.6 km s -1 . This implies 
that the FWHM, calculated as y/8 In 2a v , is of the order 
of ~ likms" 1 , a valu e close to the observed 1.2 kms -1 
by lOnishi et al.1 (|1999D for the Pipe. The discrepancy 
with the analytic result, however, is probably due in part 
to the more complex velocity field of the simulation, in- 
cluding the global collapse. In some cases the observed 
line width is inflated because of the superposition of sep- 
arate cores along the line of sight, which can be seen in 
systems with two peaks of emission (see also last column 
in Table [1]). In any case the simulations are in reasonable 
agreement with the Onishi et al. observations. 

While there are other possible mechanisms to gen- 
erate the observed velocity dispersions in molecular 
gas (one of them would be MHD turbulence, see e.g. 
iMcKee fc Zweibell (|1995| ) for a detailed discussion), we 
point out that in a cloud of many Jeans masses, of non- 
uniform density and of irregular geometry, it will be im- 
possible to avoid having significant gravitational accel- 
erations locally. These in turn will produce motions on 
the order of 1 km s _1 (see Fig. 3 and $2]), unless there 
is a special forcing or velocity field that would prevent 
this. Yet this forcing not only would have to be chosen 
so cleverly in time and space as to support just the right 
region that is about to collapse, but it also would have 
to prevent the collapse while reproducing the supersonic 
motions in the environment of the filament. Thus while 
we cannot say that there is no other source for the ob- 
served gas motions, we do argue that at least part of 
these motions are likely to be gravity-driven. 

5. DISCUSSION & CONCLUSIONS 

The results presented in this work support the 
idea that molecular clouds are in a general state of 
global collapse, su g gested more than 35 years ago by 
iGoldreich fc Kwanl (|1974h . Although such global col- 
lapse of irregular structures develops internal turbulence 
at some level, the large linewidths in MCs are caused 
mainly by the large-scale systematic inward motions in 
this scenario. Given the irregularities, angular momen- 
tum conservation and the high Reynolds numbers, it 
would actually be outright surprising if the collapse did 



n ot generate some "turbulenc e" . 

IZuckerman fc Evansl (|1974l ) suggested that such large- 
scale collapse would result in a star formation rate much 
higher than observed, mandating some mechanism of 
cloud suppo rt to render star formation inefficient (see 
lMcKed["l999l for a summary). Yet, analytical and numer- 
ical s t udies (e.g..|Burkert fc Hartmannl 2004 iField et al.l 



20081: IHeitsch et all l2008bt IHeitsch fe Hartmannl 120081: 
Vazauez- Semaden i et a T I2007L 120081 12009) demonstrate 
that MCs, while globally collapsing, are highly suscepti- 
ble to strong fragmentation. While this fragmentation is 
to some extent a consequence of a combination of ther- 
mal and dyn amical instabilities d uring the cloud forma- 
tion process (|Heitsch et al. 2008a), it is mainly due to the 
fact that the cloud geometry is finite, leading to gravita- 
tional edge focusing, i.e. non-line ar gravitational accel- 
eratio ns as a function of position ([Burkert fc Hartmannl 
l2004f) . In other words, the existence of non- linear grav- 
itational accelerations as a function of position, leading 
to a rapid piling up of material as well as local fragmen- 
tation, allows local collapse to proceed faster than global 
collapse. And it is precisely because small, high-density 
structures are rapidly developing, that the mass involved 
in the densest regions has to be small, ensuring a small 
star formation efficiency. In other words, the observed 
low star formation efficiency seems to be the testimony 
of the importance of non-linear acceleration at particular 
places while the global collapse occurs. In this connec- 
tion it is worth pointing out that the major region of star 
form ation within the Pipe (with 15 YSOs; ICovev et al.l 
|2009() . the B59 cloud, lies at one end of the filament, 
which is a preferred locus fo r gravitational edge focusing 
(jBurkert fc Hartmannll2004h . 

While there are clearly a number of uncertainties in 
both the analytical and numerical calculations and in 
the precise observational quantities, a simple model of 
gravitational infall toward the filament clearly can ac- 
count for the observed non-thermal velocity dispersions 
in the Pipe. The gravity-driven model also does pre- 
dict small (core-to- core) velocity di s persions consistent 
with observations (iHartmannl [2002; Wa lsh et al.l I2004L 
I2007t iKirk et al.ll2009D . This agreement is possible be- 
cause the "turbulence" is not being continuously driven 
by an ad hoc force but is the result of (global) gravita- 
tional acceleration. The simulation results also suggest 
that while gas is falling onto the filament supersonically, 
the velocity dispersions in the dense post-shock gas are 
subsonic, consistent with the more detailed treatment o f 
core formation in post-s hock gas bv lGomez et al.l (2007); 
IGong fc Ostrikerl (|2009f ). Yet the current simulation is 
not sufficiently resolved to support more detailed state- 
ments about t he de nse cores. 

lLada et all (|2008h suggested that the confining pres- 
sure for many Pipe cores was the result of the "weight" 
of the cloud. Here we refine this suggestion by point- 
ing out that a static pressure estimate is not strictly ap- 
propriate for a cloud in supersonic motion. The ram 
pressure that we infer is, however, another aspect of 
the same physical mechanism - gravitational accelera- 
tion. It is also worth noting that, in this interpreta- 
tion, the non-thermal velocity dispersion in the cores and 
surrounding regions are dominated by collapse motions 
and thus do not provide pressure support against grav- 
ity (|Ballesteros-Paredes et al . 1999b; Ballcster os-Paredesl 
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120061 IVazguez-Semadeni et all [20081. Nor does in this 
picture a driving source for turbulence within molecular 
clouds appear to be necess ary (see also lField et al.ll2008t 
iHeitsch fc Hartmann|[200l ). 

The numerical simulations, with an average infall ve- 
locity in the y-direction of ~ 0.4 km s -1 and a density 
of ~ 10 3 cm -3 at the filament, imply ram pressures of 
order P/k ~ 2 x 10 cm" K. Scaling this result from the 
~ 1000M o filament of the simulation to th e ~ 3000M W 
estim ated for the Pipe region from 13 CO (jOnishi et al.l 
119991) results in a predicted pressure 3 times larger, or 
about 6 x 10 4 cm -3 K. This is in reasonable agreement 
with the confining pressures for Pipe cores estimated by 
lLada et alJ (12003) of P/k ~ 8 x 10 4 cm' 3 K. 

The core mass f u nction (CMF) inf e rred for the Pipe 
(1 Alves et al.l l2007t lLada et alJ l2008t iRathborne et all 
2009) peaks at a mass well above the typically-estimated 
peak of the stellar initial mass function (IMF). In the 
gravity-driven picture, translating the CMF at an in- 
stant of time to an IMF is complicated because the 
model of gravitationally-accelerated infall implies that 
the filament - and th us the cores - are accreting mass 
(see lClark et al.l[2007l for a related argument). If we set 
the ram p ressure of in f all eq ual to confining pressure in- 
ferred bv lLada etall (|2008l ), and use v = OJIkms" 1 
from the analytic model, the implied density of infalling 
material p ~ 1.9 X 10~ 21 gcixr 3 or N(H 2 ) ~ 480cm" 3 . 
The mass accumulation rate over a cylindrical region 
of radius R = 0.1 pc and length z — 0.2 pc is then 
2irRzpv ~ 2.7 x lO -6 M yr -1 . While this estimate 
is obviously sensitive to the adopted parameters, it is 
clear that the lower-mass cores can quite plausibly dou- 
ble their mass over timescales of less than 1 Myr (see also 
iGomez et al]l20"07t iCong fc Ostrike3l20"?)flh . In fact, ob- 
servations of starless cores are starting to indicate that 
quiescent, c oherent dense star less cores do accrete ac- 
tively (e.g.. ISchnee et all l2007f ). Alternatively, the nu- 
merical simulation sho ws a doubling of mass in core s over 
a timescale of 1 Myr (|Heitsch fc Hartmannl [200 8[) , con- 
sistent with this estimate. 



The long-range nature of gravity makes it very difficult 
to avoid the generation of supersonic velocity fields and 
the production of filaments and ot her dense structures in 
clouds with many Jeans mass es (iBurkert fe Hartmann l 
20041: iBate et all I2002L 120031: IVazauez Semadeni et al.l 
20071 IHeitsch et al.ll2008bt ISmith et al.ll2009h . 

The mass, spatial extent, and typical column densities 
of the Pipe are not very different from those of the Tau- 
rus molecular cloud. It is conceivable that the Pipe is 
representative of how Taurus would have appeared 1 - 
2 Myr ago. The average filament line density in Tau- 
rus is estimated to be roughly twice the critical value 
(jHartm ann 2002) , which is consistent with the active star 
formation in the region. Blindly applying the mass infall 
rate calculated above to Taurus would imply a buildup 
from ~ 15M pc _1 to the observed ~ SOMqpc^ 1 on a 
timescale of 1 Myr. 

Finally, one may consider the future evolution of the 
Pipe nebula with regards to the star formation efficiency. 
While we have shown that besides the gas dynamics, 
the gravity-driven model can naturally explain the origin 
of the strong fragmentation as a key ingredient for the 
observed low star formation efficiency, we have not dis- 
cussed how the remaining molecular gas can be prevented 
from eventually collapsing onto the dense filament. In 
low-mass star-forming regions this could be achieved 
by e.g. magnetic support of the diffu se molecular gas 
(|Hever et all 120081: iPrice fc Batdl2009D. or bv G alactic 



tidal disruption (jBallesteros-Paredes et aT1l2009al |bl). 



We thank the referee for comments and questions that 
helped to improve our presentation. This work was sup- 
ported in part by NSF grant AST-0807305 and by the 
University of Michigan. Computations were performed 
at the National Center for Supercomputing Applications 
(AST 060034). J.P.-B. acknowledges support from grant 
UNAM-PAPIIT INI 10409. This work has made use of 
NASA's Astrophysics Data System. 



REFERENCES 



Alves, J., Lombardi, M., & Lada, C. J. 2007, A&A, 462, L17 
Andre, P., Belloche, A., Motte, F., & Peretto, N. 2007, A&A, 472, 
519 

Audit, E. & Hennebelle, P. 2005, A&A, 433, 1 
Ballesteros-Paredes, J. 2006, MNRAS, 372, 443 
Ballosteros-Paredes, J., Gomez, G. C, Loinard, L., Torres, R. M., 

& Pichardo, B. 2009a, MNRAS, 395, L81 
Ballesteros-Paredes, J., Gomez, G. C, Pichardo, B., & Vazquez- 

Semadeni, E. 2009b, MNRAS, 393, 1563 
Ballesteros-Paredes, J., Hartmann, L., & Vazquez-Semadeni, E. 

1999a, ApJ, 527, 285 
Ballesteros-Paredes, J., Klessen, R. S., Mac Low, M.-M., & 

Vazquez-Semadeni, E. 2007, in Protostars and Planets V, ed. 

B. Reipurth, D. Jewitt, & K. Keil, 63-80 
Ballesteros-Paredes, J., Vazquez-Semadeni, E., & Scalo, J. 1999b, 

ApJ, 515, 286 

Bate, M. R., Bonnell, I. A., & Bromm, V. 2002, MNRAS, 332, L65 
— . 2003, MNRAS, 339, 577 
Bernes, C. 1979, A&A, 73, 67 
Burkert, A. & Hartmann, L. 2004, ApJ, 616, 288 
Clark, P. C, Klessen, R. S., & Bonnell, I. A. 2007, MNRAS, 379, 
57 

Covey, K. R., Lada, C, Muench, A., & Roman-Zuniga, C. 2009, in 
Bulletin of the American Astronomical Society, Vol. 41, Bulletin 
of the American Astronomical Society, 399 

Dobbs, C. L. & Bonnell, I. A. 2007, MNRAS, 376, 1747 



Elmegreen, B. G. 1979, ApJ, 231, 372 
— . 2007, ApJ, 668, 1064 

Field, G. B., Blackman, E. G., & Keto, E. R. 2008, MNRAS, 385, 
181 

Forbrich, J., Lada, C. J., Muench, A. A., Alves, J., & Lombardi, 

M. 2009, ArXiv e-prints, 0908.4086 
Goldreich, P. & Kwan, J. 1974, ApJ, 189, 441 

Gomez, G. C, Vazquez-Semadeni, E., Shadmehri, M., & 

Ballesteros-Paredes, J. 2007, ApJ, 669, 1042 
Gong, H. & Ostriker, E. C. 2009, ApJ, 699, 230 
Hartmann, L. 2002, ApJ, 578, 914 

Hartmann, L., Ballesteros-Paredes, J., & Bergin, E. A. 2001, ApJ, 
562, 852 

Hartmann, L. & Burkert, A. 2007, ApJ, 654, 988 

Heitsch, F. & Hartmann, L. 2008, ApJ, 689, 290 

Heitsch, F., Hartmann, L. W., & Burkert, A. 2008a, ApJ, 683, 786 

Heitsch, F., Hartmann, L. W., Slyz, A. D., Devriendt, J. E. G., & 

Burkert, A. 2008b, ApJ, 674, 316 
Heitsch, F., Slyz, A. D., Devriendt, J. E. G., & Burkert, A. 2006a, 

MNRAS, 373, 1379 
Heitsch, F., Slyz, A. D., Devriendt, J. E. G., Hartmann, L. W., & 

Burkert, A. 2006b, ApJ, 648, 1052 
Hennebelle, P., Banerjee, R., Vazquez-Semadeni, E., Klessen, R. S., 

& Audit, E. 2008, A&A, 486, L43 
Heyer, M., Gong, H., Ostriker, E., & Brunt, C. 2008, ApJ, 680, 420 
Kim, W.-T., Ostriker, E. O, & Stone, J. M. 2003, ApJ, 599, 1157 



10 



Heitsch et al. 



Kirk, H., Johnstone, D., & Basu, S. 2009, ArXiv c-prints 
Kirk, H., Johnstone, D., & Tafalla, M. 2007, ApJ, 668, 1042 
Klessen, R. S., Ballesteros-Paredes, J., Vazqucz-Scmadcni, E., & 

Duran-Rojas, C. 2005, ApJ, 620, 786 
Koyama, H. & Inutsuka, S.-i. 2002, ApJ, 564, L97 
Lada, C. J., Mucnch, A. A., Rathborne, J., Alves, J. F., & 

Lombardi, M. 2008, ApJ, 672, 410 
Li, Z.-Y. 2001, ApJ, 556, 813 

Lombardi, M., Alves, J., & Lada, C. J. 2006, A&A, 454, 781 
Mac Low, M.-M., Klessen, R. S., Burkert, A., & Smith, M. D. 1998, 

Physical Review Letters, 80, 2754 
McKee, C. F. 1999, in NATO ASIC Proc. 540: The Origin of Stars 

and Planetary Systems, ed. C. J. Lada & N. D. Kylafis, 29 
McKee, C. F. & Zweibel, E. G. 1995, ApJ, 440, 686 
Muench, A. A., Lada, C. J., Rathborne, J. M., Alves, J. F., & 

Lombardi, M. 2007, ApJ, 671, 1820 
Myers, P. C. 1983, ApJ, 270, 105 

Onishi, T., Kawamura, A., Abe, R., Yamaguchi, N., Saito, H., 
Moriguchi, Y., Mizuno, A., Ogawa, H., & Fukui, Y. 1999, PASJ, 
51, 871 

Ostriker, J. 1964, ApJ, 140, 1056 

Padoan, P. & Nordlund, A. 1999, ApJ, 526, 279 

Patel, N. A., Goldsmith, P. F., Heyer, M. H., Snell, R. L., & Pratap, 

P. 1998, ApJ, 507, 241 
Price, D. J. & Bate, M. R. 2009, ArXiv c-prints 
Rathborne, J. M., Lada, C. J., Muench, A. A., Alves, J. F., 

Kainulainen, J., & Lombardi, M. 2009, ArXiv c-prints 



Rathborne, J. M., Simon, R., & Jackson, J. M. 2007, ApJ, 662, 
1082 

Schnec, S., Caselli, P., Goodman, A., Arce, H. G., Ballesteros- 
Paredes, J., & Kuchibhotla, K. 2007, ApJ, 671, 1839 
Schneider, S. & Elmcgreen, B. G. 1979, ApJS, 41, 87 
Smith, R. J., Clark, P. C, & Bonnell, I. A. 2009, MNRAS, 667 
Stone, J. M., Ostriker, E. C, & Gammie, C. F. 1998, ApJ, 508, 
L99 

Tobin, J. J., Hartmann, L., Furesz, G., Mateo, M., & Mcgcath, 

S. T. 2009, ApJ, 697, 1103 
van Dishoeck, E. F. & Black, J. H. 1988, ApJ, 334, 771 
Vazqucz-Semadeni, E., Gomez, G. C, Jappsen, A. K., Ballesteros- 
Paredes, J., Gonzalez, R. F., & Klessen, R. S. 2007, ApJ, 657, 
870 

Vazqucz-Semadeni, E., Gomez, G. C, Jappsen, A. K., Ballesteros- 
Paredes, J., & Klessen, R. S. 2009, ArXiv c-prints 

Vazqucz-Scmadcni, E., Gonzalez, R. F., Ballesteros-Paredes, J., 
Gazol, A., & Kim, J. 2008, MNRAS, 390, 769 

Walsh, A. J., Myers, P. C, & Burton, M. G. 2004, ApJ, 614, 194 

Walsh, A. J., Myers, P. C, Di Francesco, J., Mohanty, S., Bourkc, 
T. L., Gutermuth, R., & Wilner, D. 2007, ApJ, 655, 958 

Zuckcrman, B. & Evans, II, N. J. 1974, ApJ, 192, L149 



